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' Abstract. We study equilibrium and non-equilibrium properties of a two-level quantum dot 

' close to the singlet-triplet transition. We treat the on-site Coulomb interaction and Hund's rule 

coupling perturbatively within the Keldysh formalism. We compute the spectral functions and 
' the differential conductance of the dot. For moderate interactions our perturbative approach 

captures the Kondo effect and many of the experimentally observed properties. 

o 

1. Introduction 

Handling strongly interacting multilevel systems under non-equilibrium conditions is of crucial 
importance for understanding transport properties of molecules and correlated mesoscopic 
structures. Although the theoretical description of these systems is rather satisfactory in 
^ ' equilibrium [H El El HI El E] , the success of non-equilibrium methods is rather limited 0: In fact, 

- ' most of the available methods are unable to capture the Kondo physics or are rather uncontrolled. 

Perturbation theory in this regard is of particular importance: Although it breaks down for 
^ . strong interactions, for moderate interaction strengths it is able to capture the formation of the 

Q^ I Kondo resonance and the Hubbard peaks [8], scales very well with the number of orbitals, and 

■ it is therefore a promising candidate to combine with ab initio calculations. 
Motivated by experiments on lateral quantum dots |9] and carbon nanotubes [TU], here we 

focus on a particular parameter range of multilevel quantum dots, and study the so-called singlet- 
triplet transition, i.e., the transition from a singlet state of the dot to a Kondo-screened triplet 
5^ I state, driven by the presence of Hund's rule coupling. We show that the so-called interpolative 

perturbation theory (IPT) O [TTl [T21 [13] can be extended to include Hund's rule coupling, and 
it captures such basic features of this transition as the splitting of the Kondo resonance [9j . A 
detailed analysis of the transition shall be published elsewhere [13] ■ 

2. Model and theoretical framework 

We use the following two level Hamiltonian to investigate the transition, H = Hq + Hi^t, 

Ho = J2 ^"4a'^C«a<^ + Y ^i4a(^i<^ + Y *ai(4"<7^i<^ + ^-C') > (1) 
Hint = U/2Y,(i,a)^{j,a')^ianja' - JS"^ . 
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Figure 1. The symmetrized vertex and the first and second order self-energy diagrams 



Here c^^^ creates a conduction electron in the left or right lead with energy = C+^o (Mq = eV"a 

is the bias applied on lead a S (L, Kj) and spin o and is the creation operator of an electron of 
spin o", and energy e± on the dot level i G (+, — ). The ta.% denote the tunneling matrix elements 
between lead a and dot level i. Throughout this paper we focus on a completely symmetrical 
lateral quantum dot, and assume that one of the dot levels is even (+), while the other is odd 
under reflection (— ). As a consequence, the tunneling matrix elements have a simple structure: 
= = u+/\/2 and II- = —tR,- = V-/V^ [H IS]- The width of the dot-levels is given 
by F-i- = 27rpo\v±\'^, with po the density of states of the electrons in the leads. The coupling U 
denotes the on-site Coulomb interaction, and accounts for the charging energy of the dot, while 
J stands for the Hund's rule coupling which favors a ferromagnetic alignment of the total spin 

of the dot, S = i Y.iaa' 4a^dia. 

The simple Hamiltonian ([T]) describes a variety of physical phenomena, and captures, e. g. 
the underscreened [T3] and the 5C/(4) Kondo states [ini ISl dZl IE]- Here, however, we shall 
focus only to the regime with (n_|_ -|- n_) ~ 2, and the vicinity of the singlet-triplet transition 
induced by the competition of the Hund's rule coupling and the separation of the two levels, 
A = e+ - e_. 

We shall treat ([T]) by applying perturbation theory in U and J. However, before doing so, 
we separate the Hartree contribution by introducing counterterms, Hq = Hq + -ffcountj 
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We then use Hq to obtain the unperturbed Keldysh Green's functions, g'^^^, with i and i' dot- 
level labels and k and k' = i: the usual Keldysh labels, and treat -ffcount = ^{^i — ^ia)d\^di(j 
as a perturbation. 

We then need to do second order perturbation theory in ifmt- To treat both the Coulomb 
interaction and the Hund's rule coupling on equal footing, we merge them into a single interaction 
vertex, T^^^™*^ , as diagramatically shown in Fig. [TJ Then the second order self-energy can be 
expressed in time domain as follows, 
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Up to second order, the frequency dependent Green's functions can be obtained from the 
following Dyson's equation 



(4) 



where Ti^^^'^^^(uj) contains the self-energy part coming from the counterterm and the first order 
Hartree contribution, also shown in Fig. [H 



Eqs. dSj) and ^ give a complete perturbative description of the quantum dot, however, they 
depend parametrically on the so far unspecified levels, £±a- These are determined selfconsistently 
from Eqs. ([3]) and ^ by requiring that g and G give the same occupation numbers [2] 

nS? [e^a] = ^ / gZi-^)du; ^ — / GZ{u^)du; = nfj [^^<r\ ■ (5) 

Equilibrium results: For simplicity, here we focus on the electron-hole symmetric case, where 
e± = —1.5U lb A/2. First we consider the equilibrium spectral functions of the two levels: 

pM = -Ism Gl^iu) = {G>M - G<^{co)) , (6) 

the sum of which, pt{oj) = Yli a Pi'^i'^) •> plotted in Fig. E) Consider first the case J = 
(Fig. Eh). In the absence of hybridization, the ground state is highly degenerate for A = 0, and 
turning on a finite hybridization leads to the appearance of a Rondo resonance. Changing the 
value of the level splitting (A/{7) the two electrons are forced to stay on the lower level, e_, and 
a dip opens in the spectral function, corresponding to a singlet ground state. 

In Fig. [2b we analyze the effect of Hund's rule coupling, J. For J ^ A a triplet ground 
state is favored, and a Kondo effect develops, where the spin 5 = 1 of the dot is screened by 
the even and odd combinations of the conduction states. A clear signal of the singlet-triplet 
transition is that the zero bias dip of the spectral function closes and a Kondo peak develops 
on the triplet side. These features are surprisingly well captured by the simple perturbative 
calculation presented here. 

Out of equilibrium results: Equations ^ and ([5]) also provide a closed set of equations for 
the local Green's function in non-equilibrium. Solving them we can then compute transport 
properties of the dot using the Meir-Wingreen formula [19]. In our calculations the effective 
energies were fixed in equilibrium for a given parameter set, and we used these equilibrium 
values for all bias voltages. 

In Fig. [3] the current is plotted as a function of bias, Sp = pL — PR = eV, for two different 
parameter sets. In case oi J/U = and A/C/ = the current is suppressed and the linear 
conductance vanishes. This is due to the destructive interference between the even and odd 
channels, which both acquire a phase shift 6± = tt/2 in this Kondo regime. The other current 
plot (red crosses) corresponds to a dot with singlet ground state. The differential conductance 
displays a dip at zero bias (see inset of Fig. [3]) , which is a clear fingerprint of the non-equilibrium 
singlet-triplet transition. 
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Figure 2. Total spectral function, pt{^) = o- /Ojo-('^)- Evolution of (a) spectral functions as 
a function of level splitting. A, (b) spectral functions as a function of Hund's rule coupling, J. 
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Figure 3. Current versus dimensionless bias, 6fi/U = eV/U, for a singlet ground state 
[J/U = 0, A/U = 1.2) and for the highly degenerate ground state {J/U = 0, A/U = 0). 
Inset: differential conductance for J/U = 0, A/U = 1.2. In all cases r+ = r_ = F = 0.564 U. 



3. Conclusions 

In summary, we find that self-consistent second order perturbation theory captures qualitatively 
the transport properties of a two-level quantum dot throughout the singlet-triplet transition, 
and the splitting of the Kondo resonance. 
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